Entrainment of Coupled Oscillators on Regular Networks by Pacemakers 
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We study Kuramoto oscillators, driven by one pacemaker, on d- dimensional regular topologies with 
nearest neighbor interactions. We derive the analytical expressions for the common frequency in the 
case of phase-locked motion and for the critical frequency of the pacemaker, placed at an arbitrary 
position in the lattice, so that above the critical frequency no phase-locked motion is possible. We 
show that the mere change in topology from an open chain to a ring induces synchronization for 
a certain range of pacemaker frequencies and couplings, while keeping the other parameters fixed. 
Moreover we demonstrate numerically that the critical frequency of the pacemaker decreases as a 
power of the linear size of the lattice with an exponent equal to the dimension of the system. This 
leads in particular to the conclusion that for infinite-dimensional topologies the critical frequency 
for having entrainment decreases exponentially with the size of the system, or, more generally, with 
the depth of the network, that is the average distance of the oscillators from the pacemaker. 

PACS numbers: 05.45.Xt, 05.70.Fh 



I. INTRODUCTION 



Synchronization is an ubiquitous phenomenon, found 
in a variety of natural systems like fireflies Q], chirping 
crickets 2], or neural systems, but it is also utilized for 
artificial systems of information science in order to en- 
able a well-coordinated behavior in time Q. As an im- 
portant special case, the coordinated behavior refers to 
similar or even identical units like oscillators that are in- 
dividually characterized by their phases and amplitudes. 
A further reduction in the description was proposed by 
Kuramoto for an ensemble of oscillators 4], after Win- 
free had started with a first model of coupled oscillators 
0. Within a perturbative approach Kuramoto showed 
that for any system of weakly coupled and nearly identi- 
cal limit-cycle oscillators, the long-term dynamics is de- 
scribed by differential equations just for the phases (not 
for the amplitudes) with mutual interactions, depending 
on the phase differences in a bounded form. It is this 
model, later named after him, that nowadays plays the 
role of a paradigm for weakly and continuously interact- 
ing oscillators. In a large number of succeeding publi- 
cations the original Kuramoto model was generalized in 
various directions, for a recent review see HJ. In par- 
ticular the natural frequencies were specialized in a way 
that one oscillator plays the role of a pacemaker with fre- 
quency higher than the natural frequencies of all other 
oscillators QUI- Pacemakers play an important role for 
the formation of patterns in Belousov-Zhabotinsky sys- 
tem . Special families of wave solutions of the phases 
arise as a consequence of dynamically created pacemak- 
ers [Tjj. Moreover pacemakers are important for the 
functioning of the heart [j| and for the collective behav- 
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ior of Dictyostelium discoideum [ijj, as well as for large- 
scale ecosystems [l^. In 0] it is shown that a single 
periodic pendulum oscillator can entrain or at least dras- 
tically influence the dynamics of all chaotic pendula on 
two-dimensional lattices. Furthermore in Q the role of 
pacemakers on complex topologies was analyzed in order 
to understand the functioning of the neural network at 
the basis of the circadian rhythm in mammals. 

In this paper we consider a system of Kuramoto oscil- 
lators, coupled with their nearest neighbors on various 
regular lattice topologies, and driven by a pacemaker, 
placed at an arbitrary site of the lattice (section [nj . In 
particular, we analytically derive the common frequency 
of phase-locked motion in case of generic networks (in 
particular for d-dimensional regular lattices with open or 
periodic boundaries) in section IIIII Locked phases will 
be also called phase entrainment throughout this paper. 
We also analytically derive the upper bound on the ab- 
solute value of the ratio of the pacemaker's frequency to 
the coupling strength in case of one-dimensional regular 
lattices (section. llVjl. In section we consider higher- 
dimensional lattices and extend the results obtained for 
d = 1 to any dimension d > 2 of the lattice. We find 
that the range of pacemaker frequencies for which one 
obtains synchronization, the so-called entrainment win- 
dow, decreases with an inverse power of the linear size N 
of the lattice with an exponent given by the dimension 
d of the lattice. This leads to the conclusion that the 
entrainment window of an infinite-dimensional network 
decreases exponentially with its linear size N if the pace- 
maker is asymmetrically coupled to the other oscillators. 
This conclusion is supported by our analysis of coupled 
oscillators on a Cayley-tree, a topology that amounts to 
an infinite-dimensional regular lattice. These results con- 
firm the results recently obtained by Kori and Mikhailov 
8] for random network topologies. For random topolo- 
gies the entrainment window decays exponentially with 
the so-called depth of the network, that is the average 
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distance of all other oscillators from the pacemaker. For 
our regular topologies, the linear size of the networks for 
hypercubic lattices and the radius of the Cayley tree are 
proportional to the network depth. 



II. THE MODEL 

The system is defined on a regular network. To each 
node i, i — 0, . . . , N, we assign a limit-cycle oscillator, 
characterized by its phase ipi that follows the dynamics 



tpi = uj + S is Au + (1 + Si 
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A is the adjacency matrix of the system (A i} j - 
if the nodes i and j are connected and Aij 
erwise), it reflects the underlying topology of the net- 
work. Here only nearest neighbors are coupled. More- 
over, ki = J2j is the degree of the i-th node, it gives 
the total number of connections of this node in the net- 
work. Si j denotes the Kronecker delta ( S^j = 1 if i = j 
and Sij = if i ^ j ). The oscillator at position s rep- 
resents the pacemaker. Its natural frequency differs by 
Aw with respect to the natural frequency to of all other 
oscillators. Without loss of generality we set uj — 0, be- 
cause system is invariant under the transformation 
— > <Pi +u)t , Vi. Moreover the interaction of the pace- 
maker with the other oscillators can be linearly tuned 
by the parameter — 1 < e < 0. For e — the pace- 
maker is on the same footing as the other oscillators. For 
e = — 1 its interaction is asymmetric in the sense that 
the pacemaker influences the other oscillators, but not 
vice versa (the pacemaker acts like an external force). In 
natural systems both extreme cases as well as interme- 
diate couplings can be realized. The constant K > 
parameterizes the coupling strength. The phases of the 
i-th and j-th oscillators interact via the sine function of 
their difference, as originally proposed by Kuramoto. 



III. PHASE-LOCKED MOTION 

We consider the conditions for having phase-locked 
motion , in which the phase differences between any pair 
of oscillators remain constant over time, after an initial 
short transient time. First we calculate the frequency O, 
in common to all oscillators in the phase-locked state. 
Imposing the phase-locked condition (pi = fl , V i = 
0, . . . , N to system and using the fact that the sine 
is an odd function (see the Appendix A for details), we 
obtain 



O 



Aui 



(2) 



As long as il depends on the degree k s we see that the 
common frequency Q increases with the degree of the 



pacemaker. On the other hand, when the network has 
a homogeneous degree of connections, fcj = k , V i = 
0, . . . , N, Eq.((2J takes the form 



n = 



(l + e)N + l 



(3) 



It should be noticed that in this case the common fre- 
quency does not depend on the common degree k of the 
network. 

In terms of the original parameterization of the model 
the common frequency after synchronization is fl+to. 
For the derivation of Eq.s J5J and we made only use 
of the odd parity of the coupling function. The former 
results are still valid for any other odd coupling function 
/ (ipj — ipj) which is 27r-periodic and bounded. 



IV. ONE-DIMENSIONAL LATTICE 



s-1 



s+1 



N-l N 



Figure 1: One-dimensional lattice with N + 1 sites labelled 
by their coordinates. 



A. Linear chain 

Let us consider first the case of N+l Kuramoto oscilla- 
tors coupled along a chain (Figure^ with free boundary 
conditions at the positions and N (ki = 1 if i = 
or i = N, ki = 2 otherwise), the pacemaker located at 
position s, < s < N. In this case it is convenient 
to introduce the phase lag between nearest neighbors 
8i = ipi — (fi—\. Consider first the pacemaker placed 
at position < s < N. We start with the nearest oscil- 
lators to the right of the pacemaker, placed at position 
i = s + 1 

Q = K/2[ sin (-9 S+1 ) + sin (6 s+2 ) } 
^sm{8 s+ 2) = 2n/K + sm(9 s+1 ) ' 

Moving again to the right, the equation of the (s + 2)-th 
oscillator reads 

fl = K/2 [ sin {-e s+2 ) + sin (6 s+3 ) } 
^sin(e s+3 ) = 2n/K + sm(0 s+2 ) ■ 
=> sin (6> s+3 ) = AVL/K + sin (0 s+ i) 

Iteratively, we can write for each 1 < j < JV — s 



sin (6 s+j ) = 2 (j - 1) — + sin (6» s+1 ) 

A 

In particular when j = N — s we have 



(4) 



sm(0jv) = 2(N- 



1)_ +S in (0 a+1 ) 



3 



but also at the boundary 

Q = Ksin(-0 N ) . 

From the last two equations we can simply determine the 
value of sin (6 s +i) as function of s, N and Q. Substituting 
this value into Eq.QJ, we obtain 



sm(6 s+j ) = 2(j-l)-| + (2S-27V + 1) 



Q 
K 



(5) 



When fi > (fi < 0), Eq.© is always negative (posi- 
tive) and has its minimum (maximum) value for j = 1. 
This means that when the pacemaker succeeds in "con- 
vincing" its nearest neighbors to the right to adapt his 
frequency, all the others to the right do the same. Now, 
the absolute value of the critical threshold can be cal- 
culated by only using the fact that the sine function is 
bounded ( | sin (9) | < 1 ) and using the expression for ft 
as a function of e and AT, as it is derived in the Appendix 
A in Eq-O 



Auj 
~K 



(l + e)N-e 
27V - 2s - 1 



(6) 



Eq.@) yields the bound for oscillators to the right (R) of 
s to approach a phase-locked state. Following the same 
procedure, but moving to the left of the pacemaker, we 
find 



Auj 



(l + e)iV-£ 
2s - 1 



(7) 



as bound for oscillators to the left (L) of the pacemaker 
to synchronize in a phase-locked motion. Since we are 
interested in a state with all oscillators of the chain being 
phase-entrained, we need the stronger condition given by 



(8) 



Auj 




R 


Auj 


L 


Auj 






= min 










K 


c 




K 


C 


K 


C - 



For the pacemaker placed at the boundaries s 
s = N, using Eqs.JU and HA2jl. we obtain 



Aw _ (1 + e) (2N - 1) + 1 

~K „ ~ 2N - 1 



or 



(9) 



B. Ring topology 

If we close the chain of AT + 1 oscillators to a ring, 
fcj = 2 , Vi, the derivation of the upper bound on the 
pacemaker's frequency \Auj/K\ c proceeds in analogy to 
that of Eq.(jnj). Using Eq.JSJ), the final result is then given 
as 



Auj 
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Figure 2: Critical threshold \Auj/K\ c as function of the po- 
sition s of the pacemaker on a one- dimensional lattice. For 
further details see the text. 



In all former cases, for A^ — * oo, the critical threshold 
\Auj/K\ c goes to (1 + e) [or values proportional to 
(1 + e)]. Differently, the common frequency goes to lu 
for > e > — 1 and A^ — » oo, while f2 goes to ui + Auj for 
e = — 1 and A^ — > oo. Therefore the symmetric coupling 
of the pacemaker to the rest of the system favors the 
synchronizability of the system, while it can no longer 
synchronize for a completely asymmetric coupling of 
the pacemaker (e = —1). This is plausible as it must 
be easier for A^ oscillators to convince one pacemaker 
to follow them (case e > —1) than the opposite case, 
in which the pacemaker must convince A^ oscillators to 
follow it (case e = —1). For increasing system size, the 
latter case becomes impossible, while the former is still 
possible. 



C. Topological switch to synchronization 

The numerical results of this paper are obtained by 
integrating the set of Eq.slQJ with the Runge-Kutta 
method of fourth order (dt = 0.1). The numerical value 
of the critical threshold \Auj/K\ c is evaluated with an 
accuracy of 5 • 10~ 5 . As one knows from Q, in case of 
regular topologies, stable solutions with different winding 
numbers are possible. To avoid different winding num- 
bers, we always choose homogeneous initial conditions 
(also the distribution of phases would be different in par- 
ticular in the synchronized state). 

We summarize the results obtained so far in Figure |2] 
The analytical results for \Aui/K\ c are represented by 
lines for the open chain [dotted line from Eq.@ and 
dashed line from Eq.Q], for the ring [full line from 
Ea. (|10(l ] and by crosses [Eq.||3}] in case of an open chain 
with the pacemaker at the boundaries, while the circles 
(open chain) and squares (ring) represent numerical data 
that reproduce the analytical predictions within the nu- 
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merical accuracy. All results are obtained for N = 20 and 
e = 0, they are plotted as a function of the pacemaker's 
position s that only matters in case of the open chain. 
The horizontal line obviously refers to the ring, the two 
branches (left and right), obtained for the chain, cross 
this line when the pacemaker is placed at s = 10 in the 
middle of the chain. When the pacemaker is located at 
the boundaries s = and s = 20, we obtain two isolated 
data points close to the horizontal line. 
Let us imagine that for given N and e the absolute value 
of the pacemaker's frequency |Aw| and the coupling K 
are specified out of a range, such that the ratio is too 
large to allow for phase-locked motion on a chain, but 
small enough to allow the phase-entrainment on a ring. 
It is then the mere closure of the open chain to a ring that 
leads from non-synchronized to synchronized oscillators 
with phase-locked motion. Therefore, for a whole range 
of ratios |Au;/i<T|, no finetuning is needed to switch to 
a synchronized state, but just a simple change in topol- 
ogy, the closure of a chain to a ring. Because this closure 
may be much easier feasible in real systems than a fine 
tuning of parameters to achieve synchronization, we be- 
lieve that this mechanism is realized in natural systems 
and should be utilized in artificial ones. In our numeri- 
cal integration we simulated such a switch and plot the 
phase portrait in Figure The phases tpi as function of 
time are always projected to the interval [0, 2ir): we use a 
thick black line for the phase of the pacemaker and thin 
dark-grey lines for the other oscillators. In this concrete 
numerical simulation with T ■ dt = 4000 integration steps 
altogether, we analyzed a one-dimensional lattice of Ku- 
ramoto oscillators with N — 6, s = 2, Aoj/K — 1 and 
e = 0. In the time interval from to "ON" (T/3) we see 
a phase evolution with different slopes and varying with 
time. Moreover, the pacemaker and the left part of the 
system (oscillators i = 0, 1) have larger frequencies than 
the right part of the system (oscillators i = 3,4,5,6). 
At the instant "ON" we close the chain, passing to a 
ring topology, the system almost instantaneously reaches 
a phase-locked motion (all phases moving with the same 
and constant frequency). At time "OFF" (2T/3) we open 
the ring, again, the system then behaves similarly to the 
first phase, i.e. for t £ [0, ON). Furthermore it should be 




time 



Figure 3: Phase portrait of a system of Kuramoto oscillators 
on a one-dimensional lattice. This figure shows how it is pos- 
sible to switch "ON" and "OFF" synchronization via a simple 
topological change of the system, passing from a linear chain 
to a ring and vice versa. A more detailed description is given 
in the text. 

noticed from Figure [3 that it is also favorable to put the 



pacemaker at the boundaries of an open chain to facili- 
tate synchronization. For d = 1 the pacemaker then has 
to entrain only one rather than two nearest neighbors so 
that the range of allowed frequencies |Aw| increases. 

V. HIGHER DIMENSIONS 

All of our results obtained so far extend qualitatively 
to higher dimensions d, when system Q is placed on a 
hypercubic lattice with (Nj + 1) oscillators in each direc- 
tion j, so that we have an ensemble of Hj =1 (Nj + 1) 
Kuramoto oscillators, where the i-th oscillator's posi- 
tion is labelled by a d-dimensional vector i, with < 
ij < Nj , V j = 1, . . . , d. If the condition for having a 
phase-locked motion is satisfied, the system of oscillators 
reaches a common frequency still given by Eq.Q. Such 
condition now is satisfied for \Au/K\ < \Au>/K\ c , the 
critical ratio for the pacemaker's frequency at position 
s = (si, . . . , Sd)[see Figure 0]. In the simplest case when 

| Aco / K | c ______ f\ 

—> 1.0084 r ,X""\ ^^^S^X/f/>, 




Figure 4: Same as Figure but in d = 2 dimensions. The 
linear size of the system is N = 10. The arrow indicates the 
value of the critical threshold for periodic boundary condi- 
tions. 

the lattice has (N +1) sites in each direction, we can 
qualitatively extend all the previous results, obtained so 
far for d = 1, to any dimension d > 2. For example, 
the "closure" of the open boundaries of the lattice to a 
torus in d dimensions favors synchronization of the sys- 
tem. We checked this numerically for d = 2 and N = 10 
[Figure E). Except for the central node at s = (5,5), the 
critical threshold in case of open boundary conditions 
lies always below that for periodic boundary conditions. 
Moreover, it is natural to assume that the former re- 
sults, obtained for one-dimensional lattices, in the case 
of periodic boundary conditions extend to d-dimensional 
lattices by replacing N in Eq.JTOj! by (N + l) d - 1 
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Figure 5: Entrainment window for higher dimensional lat- 
tices. We plot the critical threshold \Au)/K\ c — (1 + e) that 
is independent on e as it is seen from Eq. llH . Fhe main plot 
shows the dependence of the entrainment window on the di- 
mension d, while the inset shows its dependence on the linear 
size N of the lattice. Numerical results are plotted as symbols. 
They fit perfectly with our predictions according to Eg . II li . 



decreases exponentially with the dimension d of the lat- 
tice and as a power of the linear size of the lattice with 
exponent d. This conjecture is supported by the numer- 
ical results as documented in Figure For this reason 
we expect that for infinite-dimensional systems (d — > oo) 
the entrainment window decreases exponentially with the 
linear size of the system if the pacemaker is asymmetri- 
cally coupled with e = — 1. 

In order to verify this conjecture, we study a system of 
limit-cycle oscillators placed on a Cay ley tree [Figure 0. 
Cayley trees have z branches to nearest neighbors at each 
node, apart from nodes in the outermost shell, where the 
number of nearest neighbors is z — 1. Cayley trees are 
infinite-dimensional objects in the sense that the surface 
of the system (i.e. the number of nodes in the outermost 
shell at distance R from the the center) is proportional 
to the volume of the system (i.e. the total number of 
nodes in the tree) . The set of evolution equations is still 
given by Eq.|T]l. with the adjacency matrix A specialized 
to the underlying Cayley tree topology. For simplicity, 
we consider only the case of the pacemaker placed at the 
center of the tree (s = 0). 

When the number of branches per nodes z is two, the 
Cayley tree becomes a linear chain, with N — 2R and 
the pacemaker placed at position s = N/2 = R. Obvi- 
ously for z = 2 we find the same results as given in the 
previous section (see the Appendix B). 
The infinite dimensionality of the tree shows up for z > 2. 
We obtain 




Figure 6: Cayley tree with z = 3 branches and radius R — 5. 



as common frequency for the phase-locked motion and 



Auj 
~K 



2(z-l) u -z 



+ (1 + 



(13) 



Q = Auj 



z-2 + (l + e) [2(z -!)*-«] 



(12) 



as critical ratio for obtaining the phase-locked motion. 
Both Eq.s and (T^J are valid for R > 2. In par- 
ticular Ea. 1)13(1 tell us that for e = —1 the entrainment 
window decreases exponentially with the radius R of the 
Cayley tree (see Figure[7J). For large z > 2 it is easily seen 
that the radius R gets proportional to the depth D of the 
network (in the limit z>2 one finds D ~ R). These re- 
sults nicely confirm those recently obtained by Kori and 
Mikhailov Q , who found that in infinite-dimensional sys- 
tems such as random networks and small- world networks 
with a high number of rewired edges, the entrainment 
window decreases exponentially with the depth of the 
network. In order to connect their results to ours, we 
note that not only the radius R of the Cayley tree, but 
also the linear size N of the hypercubic lattices are pro- 
portional to the depth of these regular networks. 



VI. SUMMARY AND CONCLUSIONS 

Entrainment of Kuramoto oscillators, coupled on reg- 
ular lattices via a pacemaker is possible for large system 
sizes (N — > oo) only in the case of symmetric (e = 0) 
couplings of the pacemaker. As the pacemaker coupling 
becomes asymmetric (—1 < e < 0), synchronization be- 
comes more and more difficult, and impossible for large 
system sizes and e = — 1. However, we are not only 
interested in the "thermodynamic" limit. For finite N 
and — 1 < e < 0, we find that for a whole range of ratios 



6 



Num. Simul. (z=3) 
Num. Simul. (z=4) 

(z-2)/[2(z-l) R 




10 



+ io- 1 



s 10 " : 
< 



10 



Figure 7: Entrainment window for Cayley trees as a function 
of the radius R. Numerical results (full dots and squares) 
refer to z = 3 and 4. The full lines correspond to plots of 
\Aui/K\ c — (1 + e) , as defined in Eg. 1131 . for different co- 
ordination numbers z, from up to down 2 = 3,4, 5, 6, 10 and 
100, respectively. 



| Aw /K | it is possible to induce synchronization by a mere 
closure of a chain to a ring. The sensitive dependence of 
synchronization on the topology in a certain range of pa- 
rameters may be exploited in artificial networks and is 
-very likely- already utilized in natural systems, in which 
a switch to a synchronized state should be easily feasi- 
ble (although we are currently not aware of a concrete 
example from biological systems). If the pacemaker is 
coupled symmetrically to the other oscillators (e = 0), 
the entrainment window stays finite in the large- N limit, 
but the common frequency approaches zero for N — > oo. 
In the other extreme case, if the pacemaker is coupled 
asymmetrically to the rest of the system (e = —1), our 
main result is that the entrainment window decreases as 
a power of the depth of the network with the dimension 
d in the exponent. Extrapolating this behavior to arbi- 
trary dimension rf, we see that one of the reasons of the 
exponentially fast "closure" of the entrainment window 
in complex network topologies 8] is their effective infinite 
dimensionality. 



From the odd parity of the sine function it follows that 
0^ = and fi, = — fl g . 

Summing QAlj l over all i except i = s we obtain 



while summing I|A1|1 over all i except i = j, where j ^ s, 
we obtain 

=> # E^s ^ = ^ + (1 + 

It should be noticed that 

= iVfc s + (iV-l)^fc 4 , 

so that we can write 



(1 + e ) k i + k s 



k s Aoj 



from which Eq.J2J is implied. 

In the case of a one-dimensional lattice with open bound- 
ary conditions, we have fco = 1 = and k 3 ■ = 2 , Vj 7^ 
0, AT, so that 



2AT - 1 , if s = or s = N 
2 (A - 1) , otherwise 



The common frequency of Eq. J2J can be written as 



n 



(i+e)(2N-i)+i > if s = or S = N 
( i +e)( A-i)+i . otherwise 



(A2) 



Appendix A: COMMON FREQUENCY 

Imposing the phase-locked condition, setting uj = 
and dividing by K/ki,V i = 0, . . . , N, system Q takes 
the form 



where 



K K 



Hi := ^sin(^ i - (pi) 



(Al) 



Appendix B: COMMON FREQUENCY AND 
CRITICAL THRESHOLD FOR CAYLEY TREES 

Consider a node in the outermost shell, i.e. at dis- 
tance R from the pacemaker, reached along the path 
sr = (si, S2, S3, . . . , sr-i, sr) from the central node. 
Here Si = 1, . . . , z , while = 1, . . . , z — 1 , Vz = 2, . . . , R 
[see Figure El label the choice of branch along the path. 
For this node we can write Eq.s QJ as 

= Ksin((pg R _ 1 -(pg R ) 

=> sin((£> ?R _i -fSn) = § 
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from which we can see that all oscillators at distance R, 
independently on the path along which they are reached 
from the center, satisfy the same equation. 
Proceeding in an analogous way as before, we find for the 
shell R - 1 



<Mi-i = ft = 

= f [E«~ii sin (ipg R - tf? R _ 1 ) + 
+ sin ((fg R _ 2 - (psn-i 

= f [ - (z - !)§ + sin (^ R -2 ~ f? R -x) 
=> sin (ip gR _ 2 - tpg^) = § [z+ (z - 1)] 
Furthermore, for the shell R — 2 we obtain 
0? R -2 = = 

f [Es~! 1= i sin(<^ s - B _ 1 -Vai,_ a ) + 
+ sin(^- B _ 3 - </? s - B _ 2 

= f [ ~ ( z ~ !)§ I z + ( z - !)] + sin (W w " 
=> sin ((^ S - R _3 - v? s - fl _ 2 ) = § [« + z(z - 1) + (z - l) s 
until we arrive for R — (t + 1) at 



from which 



sm 



SI 



For i = i? — 1 , we have 



t-i 



L 9=0 



sm (<^ - ^ ) = — 



R-2 

a^fc-lJ' + Cs-l)*- 1 

L 9=0 



but also, at the center of the Cayley tree, 

= Aw + (1 + e) § £* 1=1 sin - tpg ) = 

= Aw - (1 + e) n [z Sf^-l^ + ^-l)^ 1 " 



SI 



Aw 



R-2 

E 

9=0 



1 + (! + £) zV(z-l)« + (z-l) 



R-l 



(Bl) 

For z — 2 the Cayley tree reduces to a linear chain with 
N + 1 oscillators (iV = 2i?), open boundary conditions 
and the pacemaker placed at position s = N/2. The 
common frequency given in Ea. l|Bl() becomes the same 
as in Eq.(22J|. 

For z > 2 we can rewrite the truncated geometric series 
as 



9=0 



(B2) 



and obtain after some algebra Ea. (jl2|> . 

As we have seen, all oscillators at the same distance r 

from the pacemaker satisfy the same equation. We can 

write 



sin (ip r 





<Pr-l) = —g 



R-r-1 

z ]T ( z -iy + (z-i) R - r 

9=0 

(B3) 

suppressing the path that was followed to reach the node. 
From Ea. (|12|l it is easily seen that when SI > [SI < 0] , 
Ea. (IB3|) is always negative [positive], monotone and in- 
creasing [decreasing] and takes its minimum [maximum] 
value for r = 1. Imposing the bound of the sine function 
( |sin(<p)| < 1 ) to Eq.|B3jl, with r = 1, and inserting 
Ea. (|Bl|) we obtain 



Aw 



K 



c z 



Y. R q Zo(z-l) q + (z-l) 



R-l 



(1 



For z — 2 again we obtain Eq.sJTjJ) and J7J with s — N/2. 
For z > 2 we can again use the truncated geometric series 
of Eq.JBU and obtain Eq.(|T3J). 
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